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ABSTRACT 

Five-dimensional field theories compactified on an 5 1 /Z 2 orbifold naturally include local 
brane kinetic terms at the orbifold fixed points at the tree as well as the quantum level. 
We study the quantization of these theories before the Kaluza-Klein reduction and derive 
the relevant Ward and Slavnov-Taylor identities that result from the underlying gauge and 
Becchi-Rouet-Stora symmetries of the theory. With the help of these identities, we obtain 
a generalization of the equivalence theorem, where the known high-energy equivalence 
relation between the longitudinal Kaluza-Klein gauge modes and their respective would-be 
Goldstone bosons is extended to consistently include the energetically suppressed terms in 
the high-energy scatterings. Demanding perturbative unitarity, we compute upper limits 
on the number of the Kaluza-Klein modes. We find that these limits weakly depend on 
the size of the brane kinetic terms. 
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1 Introduction 



Field theories that realize extra compact dimensions offer new perspectives [1,2] to address 
problems associated with the breaking of the electroweak gauge symmetry of the Standard 
Model (SM). One interesting aspect of these theories is that the extra components of 
the gauge fields may acquire a vacuum expectation value through the so-called Hosotani 
mechanism [2] and so play the role of the Higgs field [3]. Another novel possibility for 
electroweak symmetry breaking emerges if the gauge symmetry is broken explicitly by 
boundary conditions, e.g. by compactifying the theory on an interval [4]. 

Higher dimensional field theories, however, are not renormalizable, in the sense that 
only a finite number of counterterms would be needed to fix the ultra-violet (UV) infinities 
to all loop orders. Therefore, higher dimensional field theories should be regarded as 
effective theories, and as such, they can play a significant role in understanding the low- 
energy properties of a more fundamental theory that could, for example, be of stringy 
origin. Their consistent formulation within the context of perturbation theory should still 
respect basic field-theoretic properties, such as gauge invariance and perturbative unitarity. 

The present study focuses on five-dimensional (5D) field theories compactified on an 
S 1 /Z 2 orbifold. Orbifolding provides a viable mechanism to introduce chiral fermions after 
the so-called Kaluza-Klein (KK) reduction of 5D theories to 4 dimensions. However, a 
consistent description of orbifold field theories requires the inclusion of brane kinetic terms 
(BKTs) which are localized at the orbifold fixed points [5,6]. These new operators are local 
and need be included as counterterms at the tree-level to renormalize the UV infinities of 
local operators of the same form that are generated at the quantum level. 

In this paper we extend our earlier approach to quantization [7] to 5D orbifold theories 
that include BKTs. An important aspect of our approach is that the higher dimensional 
theory is quantized, e.g. within the framework of generalized gauges, before the KK 
reduction. Such a gauge-fixing procedure is not trivial, since both the gauge-fixing and 
ghost sectors should contain terms localized at the orbifold fixed points. After the KK 
decomposition and integration of the extra dimension, the resulting effective 4D theory 
coincides with the one that would have been obtained if the KK gauge fields had been 
quantized mode by mode [8] in the conventional gauge. 

One of the advantages of our 5D quantization formalism is that both the usual gauge 
and the so-called Becchi-Rouet-Stora (BRS) transformations [9] take on a simple and finite 
form. In close analogy to the ordinary case of 4D Quantum Chromo dynamics (QCD), the 
invariance of the Lagrangian under these transformations gives rise to the corresponding 
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5D Ward and Slavnov-Taylor (ST) identities [10,11]. With the aid of these identities, we 
can derive a generalization of the Equivalence Theorem (ET) [12, 13] within the context of 
5D orbifold theories [14] that takes account of the presence of BKTs and the energetically 
suppressed terms in high-energy scatterings. Such a generalization is obtained by following 
a line of argument very analogous to that for the ordinary 4D case [15] and has been 
assisted by the introduction of a new formalism of functional differentiation that preserves 
the orbifold constraints on the 5D fields. 

As any perturbative framework of field theory, higher dimensional theories should 
also satisfy perturbative unitarity so as to be able to make credible predictions for physical 
observables, such as cross sections and decay widths. Requiring the validity of perturbative 
unitarity for the s-wave amplitude, we compute upper limits on the number of the KK 
modes. The restoration of high-energy unitarity differs in our models from the one taking 
place in 5D orbifold theories without BKTs [14,16]. In the presence of BKTs, unitarity is 
achieved by delicate cancellations among the entire infinite tower of the KK modes. 

We find, however, that the upper limits established from perturbative unitarity have 
a weak dependence on the actual size of the BKTs. In particular, large values of BKTs do 
not seem to screen the effect of the bulk terms. This weak dependence of the unitarity lim- 
its on the size of the BKTs may be understood from first principles as follows. High-energy 
unitarity strongly depends on the UV structure of the theory, since at high energies dis- 
tances much smaller than the compactification radius of the theory are probed. Therefore, 
only terms related to the bulk dynamics appear to be relevant when studying perturbative 
unitarity constraints. Evidently, such considerations will be of immediate phenomenolog- 
ical relevance when constraining more realistic higher-dimensional models into which the 
well-established Standard Model (SM) can be embedded. 

The structure of the paper is as follows. In Section 2, after briefly reviewing the 
ordinary QCD case, we study the quantization of 5D orbifold field theories with BKTs 
before the KK reduction. In addition, we derive the relevant Ward and ST identities 
which are a consequence of the underlying gauge and BRS symmetries of the theory. With 
the help of these identities, we derive in Section 3 a generalization of the ET, which we 
call the 5D generalized equivalence theorem (GET). The 5D GET provides a complete 
equivalence relation between the longitudinal KK gauge modes and their respective would- 
be Goldstone bosons that takes consistently into account the energetically suppressed terms 
in high-energy scatterings. Section 4 is devoted to the computation of upper limits on the 
number of the KK modes which are obtained by requiring the validity of perturbative 
unitarity. We also analyze the dependence of these unitarity limits on the actual size of the 
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BKTs. Technical aspects of our study are presented in the appendices. Finally, Section 5 
summarizes our conclusions. 



2 Ward and Slavnov— Taylor Identities 

We will first review the basic formalism of quantization and the derivation of the Ward and 
ST identities for the conventional 4D QCD case. This will help us to set up our conventions 
and obtain insight into analogous considerations concerning quantization and derivation of 
the corresponding Ward and ST identities within the context of 5D orbifold field theories 
with BKTs. 



2.1 4D QCD 

The ordinary 4D QCD Lagrangian quantized in the gauge is given by 



£qcd = -\F; u F a ^ + £ GF + £ FP , (2.1) 



where F£ v = d^A® — d v A a ^ + gj ' abc A A% is the field-strength tensor of the gluon field A a . 
Moreover, £qf and £fp denote the parts of the Lagrangian, associated with the gauge- 
fixing scheme and the Faddeev-Popov ghosts: 

£gf = -i(FH3) a = -^(^) 2 > 

C FP = c a c b = c a (5 ab d 2 - gf abc d^A c ll )c b . (2.2) 

In the absence of £qf and £fp, the remaining (classical) part of the QCD Lagrangian £qcd 
is invariant under the usual gauge transformations 

5AI = D ab b = (5 ab d^ - gf abc A c ^)6 b . (2.3) 

Hence, the gauge- invariant part of the tree-level effective action, rL4°] = — | J d 4 x F^ u F a ^ u , 
satisfies the relation r[A*] = TL4^ + SA^], from which the following master Ward identity 
for QCD is easily derived: 

d vTT- ~ 9f abc ^zrK = . (2.4) 
SA^ 5A h ^ v ' 



There is a similarity of the Ward identity ()2.4)1 with the one derived with the background 
field method (BFM) in QCD [17], which is a consequence of the gauge-invariance of the 
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effective action with respect to gauge transformations of the background fields. Therefore, 
the results derived below from the classical part of £qcd will be valid within the BFM 
framework as well. 

After functionally differentiating (|2.4jl with respect to gluon fields and subsequently 
setting all the fields to zero, we obtain Ward identities that relate tree-level n-point corre- 
lation functions to each other. With the convention that all momenta flow into the vertex, 
these Ward identities may be graphically represented as follows: 




We now turn our attention to the complete QCD Lagrangian (j2.1j) quantized in the 
covariant gauges. Although the gauge-fixing and ghost terms break the gauge invariance 
of £qcd, the complete QCD Lagrangian £qcd is still invariant under the so-called BRS 
transformations [9]: 

sA; = Dfc b , 

sc a = -^-c b c c , (2.7) 
sc a = ——^ = ^ . 

Note that the above form of BRS transformations is nilpotent after the anti-ghost equation 
of motion is imposed. The nilpotency of the BRS transformations is an important property 
that ensures unitarity of the S- matrix operator in the subspace of the physical fields [18,19]. 
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Following the standard path-integral quantization approach, we define the generating 
functional Z of the connected Green functions through the relation 



e iZ 



J DA Dc Dc exp i J d A x ( £ QC d + J^A^ + 5V + c a D a + K^sA* + M a sc a ) 

(2-8) 

In the above, J atJ- , D and D a are the sources for the gluons, ghosts and anti-ghosts, 
respectively. As usual, we have included the additional sources K a ^ and M a for the BRS 
variations sAf, and sc a that are non-linear in the fields. 

Given the invariance of the path-integral measure under BRS transformations, it is 
not difficult to show that 

Z[j a ^, D a , D a , K^, M a ] = Z[j a ^~fd f "D a , If, D a , K^+luJ^, M a -uoD~ a ] , (2.9) 

where uo is a Grassmann parameter, i.e. u 2 = 0. An immediate consequence of ()2.9|) is the 
master ST identity for QCD 

j^-^L - r/i^L + l D a d »^ = o. (2.10) 
5K a v 5M a £ SJw 



As before, a number of ST identities can be derived from ()2.10)1 . using functional differen- 
tiation techniques. 

An alternative and perhaps more practical approach to deriving ST identities is to use 
the BRS invariance of the Green function defined by means of the time-ordered operator 
formalism. To give an example, let us consider the BRS invariance of the Green function 
associated with the fields c a (x)A b v (y)A c p (z), i.e. 

s(0\Tc?(x)Al(y)A c p (z)\0) = 0. (2.11) 

Employing the BRS transformations given by ()2.7)1 . it is straightforward to obtain the 
expression 

^(0\TA;(x)A b v (y)A c p (z)\0) - %{Q\T?(x)J(y)A%z)\Q) - d;(0\Tc*(x)Al(y)c c (z)\0) 

+ gf bde (0\T^(x)c d (y)At(y)A;(z)\0) + gf cde (0\Tc a (x)A b u (y)c d (z)A;(z)\0) = 0. 

(2.12) 



The last two terms on the LHS of ()2.12j) depend on bilinear fields evaluated at the same 
space-time point y or z. These terms do not have one-particle poles for the external 
lines attached to those points, and therefore cancel on-shell. The resulting on-shell ST 
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identity for the involved one-particle irreducible Green functions may then be represented 
graphically as follows: 



c p c p 



- * ■ .....0 



p" a "^\Jj^ - a -<••■■• - r p a -* A = 0, (2.13) 



6 ;/ 



with g 2 = r 2 = 0. It is not difficult to check the validity of the above on-shell ST identity 
for the tree-level QCD Lagrangian (|2.1|) . 

For our purpose of studying high-energy unitarity in 2 — ► 2 scatterings, more useful 
are ST identities that apply to on-shell 4-point functions, e.g. to the product of fields: 
A a v {x)A h v (y)Ap(z)A%(w). Again, starting from the identity 

s(0\TF(x)[d v Al(y)][deA;(z)][d°A d a (w)]\0) = 0, (2.14) 

we arrive at the ST identity which may be given diagrammatically by 

pa p c 

kpP„q p r a (j = • ( 2 - 15 ) 

v b ad 

As we will see in the next section and Section 3, it is the 5D analog of the above on-shell ST 
identity which lies at the heart of the proof of the GET for the 5D orbifold field theories. 



2.2 5D Orbifold Theories with Brane Kinetic Terms 

We will now study a 5D Yang-Mills theory, such as 5D QCD, compactified on an S ll /Z 2 
orbifold. As has been observed in [5,6], orbifold compactification introduces BKTs for the 
5D gluon field at the orbifold fixed points y = and y = ttR. These BKTs naturally emerge 
as counterterms at the tree level, so as to cancel the UV infinities of the same form that 
are generated by radiative corrections [3]. Here and in the following, we denote Lorentz 
indices enumerating the 5 dimensions with capital Roman letters M, N etc., with M — fi, 5 
(/i — 0, 1, 2, 3), and the fifth compact dimension with y = x 5 G (—ttR, ttR]. 
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To avoid excessive complication in our analytic results, we consider 5D QCD com- 
pactified on an S ,1 /Z 2 orbifold with one BKT at y = 0. However, our discussion can very 
analogously carry over to the general case with two BKTs at y — and y = ttR. To be 
specific, the 5D QCD Lagrangian of interest reads 

A DQ CD(x,y) = + r c 5{y))F a MN F aMN + £ 5DGF + £ 5DFP , (2.16) 

where F^ N = 3mA% — OnA^ + g 5 f abc A b M A c N is the corresponding field-strength tensor 
for the 5D QCD gluon A a Ml and r c is taken to be a positive dimensionful coupling for the 
BKT at y = 0. 1 In addition, the terms £ 5 d gf and £ 5 d fp describe the 5D gauge-fixing and 
Faddeev-Popov ghost sectors to be discussed in detail below. 

In order that the 5D orbifold theory includes ordinary QCD with a massless gluon, 
it is sufficient for the 5D gluon field A a M to obey the following constraints: 

A a M (x,y) = A a M (x,y + 2nR), A a ^x,y) = A a ^(x,-y), A a 5 (x,y) = -A a 5 (x,-y). 

(2.17) 

Observe that A® (A^) is even (odd) under the Z 2 reflection, y — > —y. Although the BKT at 
y = breaks explicitly the higher- dimensional Lorentz invariance of the Lagrangian ()2.16|) 
down to 4 dimensions, the classical part of the 5D QCD Lagrangian remains invariant 
under the 5D gauge transformations 

5A a M = Dt6 b = (5 ab d M - g 5 f abc A c M )9 b . (2.18) 

For consistency, the 5D gauge parameter 9 a (x, y) has to be periodic and even under S ,1 /Z 2 : 
a (x, y) = e a {x, y + 2ttR) and 6 a {x, y) = 6 a {x, -y). 

Given the S ,1 /Z 2 parities of the 5D gluon field A a M {x, y) and the gauge parameter 
6 a (x, y), these quantities can be expanded in infinite series in terms of orthonormal functions 
f n (y) and g n {y) that are even and odd in y, respectively. More explicitly, we have 

oo 

A;(x,y) = £^ n)/4 (x)/ n (y), 

n=0 

oo 

A a 5 (x,y) = J2 A Us^)9n(y), (2.19) 

n=l 

oo 

e a (x, y ) = Y, e U x )Uv). 

n=0 

1 If r c is negative and \r c \ < 2irR, the theory will contain undesirable negative norm states [20]. 
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An important property of the orthonormal functions f n (y) and g n (y) is that the resulting 
coefficients A9,(x) and A^ 5 (x), the so-called KK modes, are mass eigenstates of the 
effective theory, after integrating out the extra dimension. For example, in the limit of a 
vanishing BKT term, r c — > 0, one recovers the standard Fourier series expansion, where 
the functions f n (y) and g n (y) are proportional to cos(ny/R) and sm(ny/R), respectively 

In Appendix El we derive the analytic forms of the orthonormal functions f n {y) and 
g n {y) for the complete 5D orbifold theory with BKTs. To properly deal with the singular 
behaviour of the BKTs, e.g. at y — 0, we introduce a regularization method to analytically 
define f n (y) and g n (y) in the interval (— e, e), where e is an infinitesimal positive constant 
which is taken to zero at the end of the calculation. In their definition interval (—ttR + 
e, 7iR + e], the e-regularized orthonormal functions are given by 



fn{y) = " X < 

V2 Sn -°7rR 



N 

9n{y) = n x < 
V2 Sn <°-KR 



cos m n y — | m n r c sin m n y , for — ttR + e < y < —e 

cosm n y , for — e < y < e (2.20) 

cos m n y + | m n r c sin m n y , for e < y < ttR + e 

sin m n y + | m n r c cos m n y , for — nR + e < y < — e 

sin m n y , for —e<y<e (2.21) 

sin m n y — ~ m n r c cos m n y , for e < y < nR + e , 



where 



N~ 2 = 1 + r c + ir 2 R 2 r 2 c m 2 n , (2.22) 

with f c = r c /(27fR) > 0. Notice that the orthonormal functions f n {y) and g n {y) are 
gauge-independent and related to each other through 

d 5 fn = -m n g n , d 5 g n = m n f n , (2.23) 

where the derivatives are understood to be piecewise defined within the respective sub- 
intervals given in ()2.20|) and (j2.21|) . Finally, the physical masses m n of the KK gauge 
modes ^4"„) M are obtained by solving numerically the transcendental equation [6]: 

= - tan (m n 7rR) . (2.24) 

In addition to the massless solution m = 0, in the limits r c — > and r c — ► oo, (|2.24|) has 
the simple analytic solutions m n >i = n/R and m n > x = (2n — 1)/2R, respectively. 

Employing Fourier-like convolution techniques developed in Appendix C.I, we can 
calculate the Feynman rules of the effective 4D KK theory. These are listed in Appendix A. 
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For example, applying the convolution techniques to the 5D gauge transformations (|2.18p . 
we find the effective gauge transformations for the KK modes 



^ A \n)ti - d^tn) ~ 9f abc ^ ° ''° ^(Qa< ^ n > l > m ' 

m,l=0 

oo 

8Af n)5 = -mjfo-gf^ 53 V^' 1 - 5 ^ e b {m) A c il)5 A n>hm , (2.25) 

m=0, i=l 

where g = g^j 1 \j2nR is a dimensionless coupling constant. In the limit r c — > 0, these 
transformations reduce to the ones stated in [7] without BKTs, with the identifications: 
A n ,i,m = S n ,i,m and A n> / i7n = 5 n) Lm- The definitions of A n ^ m and A„ >ijm are given in 
Appendix C.2. 

A technical drawback of (|2.25|) is that the number of the effective gauge-field trans- 
formations becomes infinite. Our task of deriving Ward and ST identities will be greatly 
facilitated if we introduce a formalism of functional differentiation on an S 1 jTL 2 orbifold. 
To this end, we need first to define the (^-function relevant to our S 1 fL 2 theory with BKTs. 
An appropriate construction can be achieved by exploiting the completeness relation 

oo 

S(y - y'- r c ) = ]T [ f n (y) f n (y') + g n (y) g n (y') } . (2.26) 

n=0 

In the limit r c — > 0, the ^-function goes over to the known form: 

i(Vi r c = 0) S S(y) = ± jjLg cos (f ) . (2.27) 

n=0 V 7 

For any periodic test function h(y) defined on S , it can be shown that the defining property 
of the 5-function, 

dy [1 + r c 5{y)\ h{y) 6{y; r c ) = h(0) , (2.28) 

-ttR 

holds true. In fact, the proof of ()2.28|) is based on the orthonormality of the mass eigenmode 
wavefunctions f n (y) and g n (y). 

We are now in a position to introduce a functional differentiation formalism that 
respects the S 1 /Z 2 properties of the fields A^x, y) and A^(x, y). To be specific, functional 
differentiation on S 1 /Z 2 may be defined as follows: 

8Af.{xi,y\) 1 , r 1 ^ 

oW " 5( yi -y 2 ;r c ) + 5(y 1 + y 2 ;r c ) 5^( Xl -x 2 ), 



5Ab(x 2 ,y 2 ) 2 
SA a 5 ( Xl , yi ) 1 
6Al(x 2 ,y 2 ) 2 



-5 



ab 



<%i - y 2 ; r c ) - 8{ Vl + y 2 ; r c ) ] 5 (4 \ Xl - x 2 ) . (2.29) 
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Notice that functional differentiation with respect to A b v (x2,y2) (^.5(^2,2/2)) is even (odd) 
under y 2 — > —1/2, as it should be. 

Like in the ordinary QCD case, we may now derive Ward identities that result from 
the 5D gauge invariance of the classical part of the tree- level effective action: rL4^-] = 
— J J d 4 x J^ R dy [1 + r c 5(y)} F^ IN F a MN . Making use of the newly introduced functional 
differentiation formalism, we obtain the master Ward identity for 5D QCD 

9m -^t - 9^ bC ^M = • (2.30) 

In order to translate the 5D master Ward identity into the effective 4D one, it proves useful 
to define the decompositions of the functional derivatives 8/8A a (x,y) and S/6Af(x,y) in 
terms of the orthonormal functions f n (y) and g n {y)'- 

r- OO r- 

£/„(</) " 



X 



™^y) = (2 ' 31) 

Observe that the definitions in ()2.31|) are consistent with the functional differentiation rules 
stated in ()2.29|) . Substituting (j2.31|) into ()2.30|) and integrating over y yields the effective 
master Ward identity for the compactified theory 

d 5V + m 5V = af abc V ^/2- 1 - s ^~ 5 ^ 5 ^ 

/ 5T 5T ~ 

X I JJb A (0(U A ™>^ + JJb A {l)b /X rn,n,l 
V 0/i (m)p 0/i (m)5 

By functional differentiation with respect to KK fields, we may derive Ward identities 
among different n-point correlation functions, which may graphically be represented as 
follows: 




C p C p 




(2.33) 



rz —lSnfi—Smfi—Sifi . r bd (0 9 d HP 

+ v 2 &n,m,i 9 [I '\/\/\/\/^ + / '\/\/\/\/^ 

du c p dp bv 
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(2.34) 



+ y/2 



-i-s„ 



g [^m,n,lf abd 



Aracd 
m,n,lj 



(m) p 



d ;/ 



c 

(0? / 

(n)fc V (m)p d c d 6 



(2.35) 




(2.36) 



12 




d 

(k)r / 

oo ' 

3=0 (j) k+p \ (f) q 

\ 



d c 
(fc)r / (0? / 

i A face e i A fade e 1 
+ &l, n ,jj / + &k,n,jj / J 

(j) k+q L (m) p (j) k + r <-i (m) p 

*• bv *- bv 

(2.37) 

At this point, we should remark that A^ 5 satisfies Ward identities very analogous to those 
found for the would-be Goldstone bosons in spontaneously broken gauge theories [21]. 
Moreover, we should stress again that the above 5 Ward identities ()2.33J) - (j2.37J) will also 
hold in the case of a BFM formulation of 5D Yang-Mills theories [22]. 

We will now discuss the part of the 5D Lagrangian in (|2.16|) . associated with the 
gauge-fixing scheme and the ghost sector. Within the framework of the generalized 
gauges, all mixing terms between the KK gauge fields A? s and their would-be Goldstone 
bosons ^4" n ) 5 are absent. To incorporate this property, we proceed very analogously to [7] 
and choose the 5D gauge-fixing functional 

F[A a M ] = d»Al - id,Al. (2.38) 

However, we multiply the 5D gauge- fixing and the Faddeev-Popov ghost terms with the 
expression [l+r c 5(y)] which includes localized interactions of the same form as those present 
in the gauge-kinetic part of the Lagrangian, i.e. 

Adgf(x,z/) = - [l+r c S(y)} 1 {F[A a M }) 2 = - [l + r c 5(y)} 1 (flMJ - £Mg 

CbdMw) = [1 + r c 6(y)] c a c b (2.39) 

= [1 + r c 6{y)] c a [ 5 ah (d 2 - £ dl) - g 5 f abc ( 0M° - £d 5 A c s ^ 1 ~ 6 

where c a (x,y) and c a (x,y) are the 5D ghosts that are even under Z2. 
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An important constraint in our formulation is the condition (|B.1|) of vanishing of the 
surface terms, i.e. the fact that total derivatives of fields or product of fields with respect 
to y do not contribute to the action. Unlike the naive Fourier modes, the mass eigenmode 
wavef unctions f n (y) and g n {y) comply with this crucial constraint and hence become the 
appropriate basis of decomposition of the 5D fields [cf. (j2.19j) ]. Detailed discussion is given 
in Appendix B, where we show explicitly that upon compactification, £ 5D GF and £ 5 d fp 
lead indeed to an effective 4D theory quantized in the conventional gauge. 

It is not difficult to verify that the complete 5D QCD Lagrangian (|2.16|) is invariant 
under the 5D BRS transformations 

oA a nab b 

b M — U M c ' 
„ tabc 

sc a = - — — cV, (2.40) 



s c 



2 

FWm\ 



Following the Fourier-like convolution method developed in Appendix C.I, we derive the 
BRS transformations for the KK modes of the effective theory 

oo 

q4° - F) r a nf abc V 1 /o ~ 1 ~ <5n > ~ <Sm <°~ <5 '' A Ac b 



m,l=0 
oo 



sA a (n)5 = -m n cf n) - gf abc £ ^""''"A^^c 6 



(0 ' 



m,l=0 

cabc 00 



Sti 



gf aC [K -1-^,0-^,0-^,0 A b c 

(n) — ^~ 2^ V ^n,m,l c (m) C (i) > 

m,l=0 

1 PA[ n)li - m n A? n)5 . (2.41) 



SC (n) 



In close analogy to our discussion for the ordinary QCD case in the previous section, 
we define the generating functional Z of the connected Green functions through the relation: 

dy Adqcd (2.42) 



e lZ = I DA M DcDc exp<( i / (fx 

-ttR 



tR 



+ f dy[l + r c 8(y)} (J aM A a M + 5°C° + c a D a + K aM sA a M + M a sc a ) X . 

J-nR J J 

The response of the generating functional Z under the 5D BRS transformations gives rise 
to the following master ST identity for the 5D theory: 

jaMjZ _ a bZ^_ l D a dtl 6Z_ _ SZ = 

5K aM 5M a ^ £ 5J a^ ^ §J a5 V J 
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From ([2.43)1 . we obtain the effective master ST identity 



0JX (n) (n) 01V1 {n) S 0J ( n ) 0J (n) 



n=0 



. (2.44) 



As in the ordinary QCD case, it proves more practical to derive ST identities from the 
BRS invariance of Green functions. For instance, we may translate the obvious identity, 



s^lT^^A^.^^H^^IO) = 0, 
into the corresponding on-shell ST identity, which is graphically given by 



(2.45) 



Pi 



m 



n (n) P2 




(n) fei 



a v 




(2.46) 



Likewise, the BRS invariance of the Green function pertinent to the product of fields 

S (0|Tcf n) (x)4 n)5 (z)^ n)p HAf n)(r (n)|0) = 0, (2.47) 
gives rise to following on-shell ST identity: 



Pi 




(n) fei 



(n) P2 



/ (n) fe 2 



+ — 




+ 



(2.48) 
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As we will see in the next section, the above on-shell ST identities are important to prove 
the GET in 2 — > 2 scatterings. In particular, by virtue of these identities, we see again 
that the KK modes behave very analogous to the would-be Goldstone bosons of non- 
linearly realized spontaneously broken gauge theories [23,24]. A detailed list of Feynman 
rules, including would-be Goldstone-boson and ghost interactions, is given in Appendix A. 



3 The 5D Generalized Equivalence Theorem 

In the high-energy limit of a spontaneously broken gauge theory, the amplitude for emission 
or absorption of longitudinal massive vector bosons equals, up to a phase and possibly up to 
a renormalization scheme dependent constant, the amplitude for emission or absorption of 
the associated unphysical would-be Goldstone modes. This high-energy equivalence relation 
between the longitudinal gauge bosons and the would-be Goldstone bosons constitutes the 
famous Equivalence Theorem (ET) [12,13]. 

In the previous section, we have shown that the on-shell ST identities, which are 
required for a rigorous proof of the ET [15], are very analogous to those of a spontaneously 
broken gauge theory, if the scalar gauge modes ^4" n ) 5 are identified with the would-be 
Goldstone bosons. Such an identification has proven very important in order to establish 
the validity of the ET in 5D orbifold Yang-Mills theories without BKTs [14]. Here, this is 
extended to 5D orbifold theories with BKTs. In detail, the ET reads: 

T( A ai A ak <? — > A bl A bl <3'\ — 

where m ni are the KK masses, E is the centre of mass system (c.m.s.) energy of the 
high-energy scattering process. The parameter C is a constant that generally depends on 
the renormalization scheme [25]. However, it is C = 1 at the tree level and in certain 
renormalization schemes that maintain the Ward identities of the classical action [21,26]. 
Finally, S, S' collectively denote all particles in the initial and final state which have no 
longitudinal modes. In order for the ET ()3.1|) to hold, the underlying gauge structure of 
the theory has to ensure cancellations of all terms that grow with energy in the amplitude 
on the LHS of (|3.1|) . As we will see below, in higher-dimensional models with BKTs, such 
cancellations are much more intricate, because the complete infinite tower of the KK modes 
is involved. 

The ET is a direct consequence of the so-called generalized equivalence theorem 
(GET) [15]. The GET goes beyond (|3.1)l . by providing an exact relation between the 
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relevant amplitudes which includes the energetically subleading terms. To state the GET, 
we first note that the longitudinal polarization vector can be written as 

e£(Jfe) = k»/m n + a»{k) , (3.2) 

where a^{k) = 0(m n /E) is the remainder of the polarization vector that vanishes in the 
high-energy limit. As an example, we may consider the GET applied to the production of 
two longitudinal KK gauge modes: 

^{^ a {n)T^ L \n)T ~^ ^(n)L^-(n)x) = 

- ^(^(n)T4«)T^^(n)5^(ri)5) _ ^(^(ra)T^(n)T — * ^(n)5 a (ri)) 

— iT(A a {n - )T A b ijl - )T — > d( n )^ n ) 5 ) + T(A a ^ n - )T A h ^ n - )T — > a^a^) . (3.3) 

In (13. 3|) . a? n N indicates an amplitude, where the longitudinal polarization vector of the 
corresponding gauge boson has been replaced by its remainder vector a M . Hence, the 
corresponding amplitudes are suppressed in the high-energy limit and the ET (j3.1j) is 
recovered. 

It is straightforward to apply the GET to processes with an arbitrary number of 
longitudinal gauge modes. Specifically, the RHS of ()3.3|) may be found by performing the 
following four operations: 

(i) Write down all amplitudes that result from replacing any number of longitudinal 
vector bosons by the respective would-be Goldstone bosons. 

(ii) Replace in each amplitude the remaining longitudinal polarization vectors by a? n y 

(iii) Multiply each of the resulting amplitudes with the factors i k (—i) 1 , where k is the 
number of the would-be Goldstone bosons in the initial state and I the corresponding 
one in the final state. 

(iv) Sum the complete set of the so-generated amplitudes. 

Observe that our GET relation in ()3.3|) is consistent with the above rules. 

Let us present an explicit proof of the GET for the simple example of (j3.3J) . based on 
the ST identities discussed in the last section. Using the decomposition, e\ = —k\jm n + 
cti, with ki being the incoming momentum, the ST identity (|2.4fij) and the transversality 
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condition e • k = 0, we find diagrammatically 




tb £2 £b £2 £b <52 



where the polarization vectors multiplying the amplitudes have been explicitly displayed 
at the corresponding external leg. 

Our next step consists in making use of e 2 = —k2/m n + a 2 , the ST identity (|2.48j) . 
and the fact that k\ • a± = m n . Substituting all the above into the RHS of (|3.4jl . we obtain 




e b e fa e b 



(3.5) 

It is obvious that ghost contributions cancel, and hence ()3.5|) is nothing than (j3.3|) in 
diagrammatic form. Following the same reasoning, one can prove the GET for any other 
process. At the tree level, a proof based on the Ward identities fJ2.33|) - (j2.37|) is possible 
but tedious. 

Having established the GET and the resulting ET, we now show how the ET is realized 
in a specific tree-level calculation. Consider the elastic scattering of two longitudinally 
polarized gauge bosons, A?-. L A b ,, L — > ^( n )L^(n)L- m terms of Feynman diagrams, the 
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tree-level amplitude is given by the infinite set of diagrams 



"L . (n)pi r— 1 

!)p2_p-^L. 



(n)fci 




A- crossings 



j=0 



iT 4 + iT s + %T % + iT u 



(3.6) 



where s, t and u refer to the pertinent channels. Using the Feynman rules of Appendix A 
and the relation 

(3.7) 



2m n (3 



[(l+/3 2 K-(l-/3>£] 



with (3 = yl — 4m 2 /s, it is straightforward to calculate the individual diagrams. More 
explicitly, we find 

,2 



iTa 



A 



n,n,n,n „ a 

8mt 



u 



f abe f cde s(t - u) + f ace f bde t[ s-—)+ f ade f bce u[ s - — 



4 



t 



A 



■ T s _ o-^.o A 2 A n 2 fabe rcde S ( U l ) 



n 

u — s 
2t 

t - s 



1 + 



2m 



2\ 2 



2 ' 

s — m j 



(3.8) 



1 + 



1 + 



2m 2 /^ 



+ 



u 



t - 2u 

m? n (3 2 

u-2t 

m 2 n (3 2 



t — rrij 



u 



2 ' 

u — 



;i-c e )^, u = { Pl + k 2 f = -{l + c e )^- 



(3.9) 



2m V" 2m 2 /3 2 
where the Mandelstam variables are given by 

S=(P1+P2) 2 , t = { Pl + k 1 ) 2 

and we have used the notation eg = cos 9. 

Our task of verifying the ET simplifies significantly if the infinite sums are expanded 
in terms of s. To be precise, using 

s 



2 4 4 2 

m\ rm m, mi 
= 1 A -\ A — 

2 2 9 2 ' 

s — m, s s z s z s — mi 



(3.10) 



the infinite sum T s can be split into four different sums: T s = Tq + T| + + T 6 S , where 
the subscript indicates the power of rrij in the numerator of ()3.1()j) . Note that each sum 
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is convergent. As can be shown by an explicit calculation (see end of this section for a 
similar case), Tg includes terms of order m n /y/s, which are energetically suppressed and 
do not contribute to the leading part of the amplitude. Likewise, we can make analogous 
expansions for the t- and it-exchange graphs: 

,2 l 



t — m 2 



u 



u — m. 



2 1 
- + 



5/7)7 



1 - Co s 



Cg S 



2 1 8m 1 , , 
: + z — -5)™% + 

1 + c e s 1 + c e s 2 1 J 



1 



[l + c e 



\2 c 2 



(3.11) 



m- + 0(m b r s 



6„-3\ 



The different infinite sums in T S '*' U , T 2 s ' t,u , and T± t,u can be calculated by means of (|D.6|) 
and ()D.7|) . which have been derived by employing complex analysis techniques developed 
in Appendix B of [27]. In this way, we obtain 



jabe jcde | ^0 „2 



_Cg_ 

2m 2 n ' 



+ r e f 
+ r de f 



[c e + 3)(c e - 1) 2 1 - °o 
-s + — — — s 



32m* 



Am 2 n 



32m4 



Ami 



(3.12) 



7T S 



^2 ya&eycde 



A 



n,n,n,n 



C 6 s 2 _ 

\mi 



eg s 5ce 
6mS 4 



2 



+ 0(m n /y/a) 



(3.13) 



2 race rdfte 



^ / ace / 



A 



c e + 3){l-ce)_ 2 llcg - 3 8cg-5c* + 9 



™ * 32m 4 

3 + c e 



+ 



A. 



2(1 - ce) 1 



12m 2 n 
+ 0{m n /yf?) , 



6(1 - c e ) 



(3.14) 



iT u = ig 2 f ade f ' hce 



A 



(3-c e )(l + c ) s 2 _ llcg + 3 ^ 8^ + 5^ + 9 



" " " i; 1 32//), 1 , 

3-c 



2(1 + c e ) 



12m 2 , 
+ C(m n /v^), 



6(1 + 0,) 



(3.15) 



where 



A„ 



8N^7r 2 R 2 f 3 c m 2 n 



(3.16) 



Collecting all contributions (|3.12|) - (|3.15J) . we find that the ^-contributions in the s-, t- and 
w-channel graphs cancel against terms in T4. Terms linear in s are identical for each colour 
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factor. Thus, they vanish due to the Jacobi identity. The final result is then given by 
iT 4 + iT s + iT l + iT u (3.17) 



^9 (An, n, n,n X n 



ace rdbe 



2(c e - 1) 



r 2 4- 3 

_j_ jade jbce » ' 



2(c e + I] 

In (|3.17p . we only exhibit the leading contribution 0(1), which will be essential for checking 
the validity of the ET. 

Let us now consider the RHS of the ET (|3.1|) . The scalar scattering amplitude is 
given by 



S s («)pi / ' 

N (n)fci 

(«)P2 / 
/ 



\ («)pi x 

E n (i) ✓ 
y\/\/\/>< + crossings 

j=0 (n)p 2 , 7 \ 



(n)fc2 



(3.18) 



£ [ <I w + <7 y) + <I w 



3=0 



Using the Feynman rules of Appendix A, we obtain for the individual s-, t- and w-channel 
graphs 



iT\ 3) = ig 2 2""^''° A^ n f ace f dbe 
%Tfa = ig 2 2-^Kl^r de f ce 
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- t 


2(s 


- m 2 ) 
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— u 
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- m 2 ) 
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— s 


2(u 


- m 2 ) 



(3.19) 



Considering the analytic form of the couplings A n j jn in (|3.19)1 . it turns out that only a 
single expansion in terms of 1/s, i.e. 



s — m- 



1 + 



s — m 



2 ■ 



(3.20) 



3 3 

will be sufficient to evaluate the leading part of the amplitude. Using (jD.7j) . we can now sum 
the leading part of (|3.19|) . Up to energetically subleading terms 0(m n /y/s), the amplitude 
(I3.18J) equals (|3.17|) . This completes our check of the ET (|3.1|) . For vanishing r c , i.e. r c — > 0, 
we find £± n ,n,n,n —* 3 and X n — > 0, and so recover the result in [14] for the validity of the 
ET in 5D orbifold theories without BKTs. 
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We conclude this section by showing that the infinite sum involving the second term 
in ()3.2U|) is indeed subleading of order m n J y/s. We modify the second equation of (|D.7|) by 
taking the additional factor m 2 /(s — m3) into account. With this modification, an explicit 
calculation gives 



oo 

E 

j=0 



m~ 



n,j,n 



A 



Am 2 



n,n,n,n 



s — Am 2 ^ 



+ 327r 2 R 2 mlr 3 N"; ( 1 
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'R 2 mirt 



mi 
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s — Ami (s — Ami) 2 



(3.21) 



+ 



d,^R A miflNt s 



n c n 



Am 2 s 2 + Am^s " 4 d L4-5a[: 



- Srr^mlrX 



s + ^tanvr^v^ (s - Am 2 n ) 2 

The first two terms are obviously 0(m^/s). The third term has got poles at the spectrum 
and we therefore take the high-energy limit in a discrete manner: yfs^ = (n — \)/R, with 
n — > oo. The third term approaches the negative of the fourth term. Consequently, all 
0(1) terms cancel and (|3.21|) does contribute only subleading terms of order m n j\fs. 



4 High Energy Unitarity Bounds 

In higher-dimensional field theories, the coupling constants are dimensionful parameters. 
For example, for the case of 5D iV c -colour QCD, the coupling g 5 has the energy dimensions 
E~ l l 2 . On grounds of naive dimensional analysis, the s-wave amplitude a of a 2 — ► 2 
scattering involving 5D gluons behaves as 

a (E 5 ) ~ N c g*E s , (4.1) 

where E 5 is the c.m.s. energy in 5D. This means that at energy scales E 5 > l/(g 2 N c ), the 
transition amplitude ao exceeds 1, thereby invalidating perturbation theory. In addition, 
the fact that the asymptotic high-energy behaviour of a does not approach a constant, but 
grows linearly with the energy, is in gross violation with the Froissart-Martin bound [28]. 

In a theory compactified to 4D, the 3- and 4-gluon couplings are dimensionless, while 
the KK spectrum is infinite. In this case, the aforementioned energy upper limit due to 
unitarity violation is translated into a corresponding upper limit on the number of the KK 
modes. To be specific, one gets the very approximate upper bound [14] 

No , 1 



< 



R - g 2 N c 



(4.2) 
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Obviously, this upper bound depends non-trivially on the compactification radius R, the 
5D gauge-coupling constant g 5 and the number of colours N c . 

We will now compute explicitly the perturbative unitarity limits on our 5D orbifold 
theory with one BKT. For this purpose, let us first describe our approach to perturbative 
unitarity. As usual, we start by stating the optical theorem for a general i — > f transition 

2 1mT fi = J2 Tfj Tji ? (4-3) 

3 

where the sum over j is understood to include phase-space correction factors for each of 
the individual states in j. Decomposing the transition amplitude Tfi in terms of Legendre 
polynomials [29], 

oo 

T fi (s, c e ) = 16vr ]T (21 + 1) Pfa) a t (s) , (4.4) 

1=0 

with ai(s) = l/(327r) dcoTfiPi(co), and substituting the resulting expression (|4.4|) 
into (|4.3|) . we obtain a unitarity relation for the s-wave (I = 0) amplitudes 

Im[a ]/i = °3 ( a o]/i [ a o)*ji ■ ( 4 - 5 ) 

3 

Here, we denote with Uj = A(s, mf, mty/s the phase-space factors for 2-particle states, with 
X(x, y, z) = [x 2 + y 2 + z 2 — 2(xy + yz + zx)} 1 / 2 . 

Our approach to perturbative unitarity consists in absorbing the phase-space factors 
(jj into the definition of the transition amplitudes [ao]y, i.e. 

[a ]ij = ^(Tidj [a ]ij • (4.6) 

This enables us to cast (J4.5|) in the simple matrix form 

Im a = a Oq • (4.7) 

Since our unitarity limits will be deduced from 2 — ► 2 scatterings that involve longitudinal 
KK gauge modes, clq will turn out to be a square symmetric matrix. Therefore, a can be 
diagonalized by means of a real orthogonal matrix R as 

R T a R = a , (4.8) 

where do is a diagonal transition amplitude matrix. The unitarity relation ()4.7|) implies that 
not only oq, but also its imaginary part Imao can be diagonalized by the same orthogonal 
matrix R. Consequently, the following unitarity relation may be established: 

Ima = (Rea ) 2 + (lma ) 2 . (4.9) 
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Equation (|4.9|) directly implies the inequality 

lma < 1 , (4.10) 

which is equivalent, by means of (|4.7p . to the requirement that the largest eigenvalue of a or 
equivalently of ao (in absolute value terms) should smaller than 1. Notice that our approach 
to perturbative unitarity differs from the one applied originally for the SM case [13], in the 
fact that phase-space effects were not considered in the latter. However, such effects are 
non-negligible for the heavier KK modes, and should therefore be included consistently. 

Since the 5D transition amplitudes diverge with increasing energy [cf. ([4.1)1 ] . the 
unitarity constraint ()4.10|) will explicitly depend on the c.m.s. energy y^. For this reason, 
perturbative unitarity limits will constrain not only the maximum allowed number of KK 
modes, e.g. iV max , but also the c.m.s. energy, i.e. v^max = 4r«7v max , above which unitarity 
is violated. To avoid complications arising mainly from IR singularities in the KK gluon 
scatterings, we follow [14] and perform a coupled channel analysis (CCA), based on the 
inelastic colour- singlet processes A^ 5 A^ n ^ — »■ A \ m )^ A \ m )^-, with n ^ m and n,m < iV max 
and s = 4m^ max . Strictly speaking, our CCA relies on the physical processes A^ L A^ L — > 
^(m)L^(m)i' Du ^ we nave made use of the ET to replace the longitudinal KK gauge bosons 
A t n ) L w ith their respective would-be Goldstone bosons A?, 5 . Such a consideration is a good 
approximation for the lighter KK modes, with masses much smaller than the c.m.s energy 
a/s. For the heavier KK modes, with masses m(„) ~ i/i, our simplification may be justified 
by the fact that scatterings of those heavy states will be relatively suppressed by phase- 
space factors that occur in the transition amplitudes [a ] m defined in ()4.6|) . Nevertheless, 
our derived upper bounds should be regarded to be slightly conservative, in the sense that 
we have not taken into account 2^2 scatterings where all asymptotic KK gluon states 
have different masses. However, initial estimates convinced us that these contributions to 
the CCA are really negligible. 2 

The phase-space-corrected s-wave amplitudes [a ] n m for the 2^2 inelastic scatter- 
ing^ A (n)5 A (n)5 Ab (m)5 A tm)V &Te § iVeI1 ^ 

s 00 

[ao]nm = [ fl o']™ ' ( 4 - n ) 

where S n ^ m = \/{s — 4m 2 ) (s — ^m"^) are phase-space correction factors and [aQ^] nm are the 

2 These estimates will be further consolidated by our analysis of unitarity bounds from the decay widths 
of the KK states, which is given at the end of the section. 
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individual s-wave amplitudes mediated by the exchange of the j KK gauge mode A ,- 
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It is interesting to observe that in the high-energy limit s ^> (m Ti 
amplitudes [ajj ]nm diverge logarithmically, 
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the s-wave 
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(4.13) 



In the limit r c — > 0, the terms j = n + m and j = |n — m| dominate the sum ()4.11|) and 
the result [14] of a 5D orbifold theory without BKTs is recovered 



g 2 N c 

32tt 



2 In 



(4.14) 



In 2 //? 2 -m 2 /R 2 \ 

Because of the presence of the BKT, our CCA involves the infinite sum (14. II}) . thus making 
it technically more challenging than the corresponding analysis without BKTs. We perform 
this infinite sum and find the maximum eigenvalue A max of ao numerically. 

In Figure ^ we display contour plots in the plane (r c ,g 2 N c ), for different fixed values 
of the maximum KK-mode number iV max , which is obtained by requiring the perturbative 
unitarity constraint: |A max | < 1. The area that lies above the contour lines is excluded by 
perturbative unitarity. For the case without BKTs, f c = r c /(2irR) = 0, our upper limits 
are larger by a about of factor 2 than those presented in [14]. The origin of this difference 
may be traced to the effect of the phase-space factors a, , which have not been included in 
the analysis of [14] . In the presence of BKTs, we find that the unitarity bounds show only 
a weak dependence on the size f c of the BKTs. In particular, the maximum allowed value 
of g 2 N c thanks to perturbative unitarity increases only by about 5%. This behaviour may 
be understood by the fact that high-energy unitarity probes distances much smaller than 
the compactification radius R where the effect of the BKTs becomes less significant. 

An interesting consequence of the BKT is that the KK gauge modes are no longer 
stable. One would have naively expected that the larger the size f c of the BKT is, the 
larger the decay width of the nth KK gauge boson becomes, thereby leading to a potential 
violation of perturbative unitarity. However, we will show that this is not the case, in 
agreement with our results from the 2^2 scatterings. Perturbative unitarity in the 
decays of the KK gauge bosons gives rise to the constraint, 

1 



T n < m r 



(4.15) 
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Figure 1: Contour plots in the plane (f c = r c /{2irR),g 2 N c ), for fixed values of iV max . The 
area above the contour lines is excluded by perturbative unitarity. 



where T n is the total decay width of the nth KK gauge boson A?,. At the lowest order 
of perturbation theory, only two-body decays of A? n s into lighter KK gauge bosons will 
contribute. The total decay width T n is then the sum of the partial widths and may 
conveniently be expressed as 

1 



V X(ml,ml,rnf) \T {nA i)\ 2 , (4.16) 

n k,l=0 (k+l<n) 

where T^ n;k> ^ is the transition amplitude for decay A? n ^ — > A^^A^^, with n > k + I. 
Squaring T( n ,k,i) and averaging over initial states, we find 

\T {n , k ,i)\ 2 = ^NjpR^NnNtNi A n , M 

m n + m k + m t + 10( m n m fc + m k m 'l + m f m n ) • (4-17) 

Here, we should note that not all kinematically allowed decay channels of A?^ contribute 
to r„. For example, the transition amplitudes iTL^d, with < k < n, vanish identically. 

In Fig. 12 we present numerical estimates of the quantity T n /(2m n ), as functions of 
the BKT coupling f c and the KK number n of the decaying gauge boson Af-> . After its 
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Figure 2: Numerical values for T n /2m n as functions of r c and n 



rapid increase, the value of T n / (2m n ) saturates to values much smaller than 1 at large f c 
and/or large n. Consequently, we find as before that BKTs do not lead by themselves to 
a violation of perturbative unitarity in the decays of the KK gauge bosons. 
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5 Conclusions 



We have studied the quantization and high-energy unitarity of 5D field theories compact- 
ified on an S ,1 /Z 2 orbifold that include localized gauge- kinetic terms at the orbifold fixed 
points. These localized interactions, the so-called BKTs, emerge naturally in these theories 
as counterterms at tree-level to absorb the UV infinities that are generated at the quantum 
level from operators of the same form. 

Extending our approach to quantization in [7], we have developed a functional differ- 
entiation formalism to quantize the theory within the framework of generalized gauges, 
before the KK reduction. With the aid of this formalism, we were able to derive the master 
Ward and ST identities which is a reflection of the underlying gauge and BRS symmetries 
of the theory. The Ward and ST identities relate the divergence of a given Green func- 
tion to other Green functions of the theory. By virtue of such identities, we have stated a 
generalized form of the ET, which we call the 5D GET. The 5D GET extends the equiv- 
alence relation between the longitudinal KK gauge modes and their respective would-be 
Goldstone bosons to consistently include the energetically suppressed terms in high-energy 
scatterings. 

An important property of any perturbative predictive framework of field theory is 
perturbative unitarity. Requiring perturbative unitarity, we have deduced upper limits on 
the number of the KK modes as functions of the size of the gauge coupling and the colour of 
the 5D Yang-Mills theory. Our approach to perturbative unitarity takes into consideration 
phase space corrections due to heavier KK modes and so improves upon earlier studies on 
the same topic [14]. Because of this, the derived upper limits have been found to be weaker 
by a factor ~ 2 than those obtained in these studies. We have investigated the impact of the 
BKTs on these unitarity limits. We have shown that the unitarity limits weakly depend 
on the size of the BKTs. Such a behaviour may be attributed to the fact that at high 
energies distances smaller than the compactification radius are probed. As a consequence, 
high-energy unitarity strongly depends on the bulk dynamics of the theory, and therefore 
high-energy unitarity limits show a weak dependence on the size of the BKTs. 

Even though large BKTs may weakly affect the unitarity limits on the theory, they 
can, however, make the higher KK modes decouple from the dynamics of the lowest lying 
KK states [6], e.g. that of the gluons. Their presence may help to considerably alleviate 
the severe phenomenological constraints that apply to those theories from a full fledged 
analysis of the electroweak data [30]. It is therefore of great phenomenological interest to 
analyze in detail the prospects of probing such BKT-dominated theories at the LHC and 
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at future e + e linear colliders. 
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A Feynman Rules 



The two-point functions for the KK mass eigenstates, deduced from the classical part of 
the 5D Lagrangian (|2.1fij) . i.e. without the terms £50 gf and £50 fp, are given by: 



(n) p 
a fi b v 

in) p 
a fi b 

(n) p 



r£(P<»)) = i5 ah [- 9t ,4p 2 -ml) +p, Pu ] 

5 ab m n p^ , 



(A.l) 

(A.2) 
(A.3) 



The two-point function (jA.2j) refers to the non-vanishing mixing between vector and scalar 
modes that occurs in the absence of the gauge- fixing term £ 5D GF discussed in Section 12.21 
In the presence of £sdgf ; however, the two-point functions take on the form: 




r£(PW. ?(»)) = ti ab [ - 9,u {P 2 - m 2 n ) + (1 - i) ^ ] , (A.4) 



Ifc(P(»), <?(»)) = "^V - 



Correspondingly, the propagators for the vector, scalar and ghost KK modes read: 



(n) p 



(A.5) 
(A.6) 



D£(p { n)) 



i8 ab 


P 2 


— m^+ie 




i5 ab 




- t,m 2 n + ^ 




-i5 ab 


P 2 


- t™ 2 n + ie 



g^v + (1 - 



PiiPv 



p 2 - £m 2 n 



(A.7) 
(A.8) 
(A.9) 



Using the definitions introduced in Appendix Of or the effective gauge- coupling coefficients 
A nm> i, &. n ,m,i etc., the Feynman rules for the interactions of the KK mass eigenstates are 
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given by 



a p 




c p 



1 — <5n,0 — 5m, — Sl.o 

[g^( k - p)p + 9pM ~ k )u + 9u P {p - q)p] 



b ft 
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\ (m) p 
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l-l ! 
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a 



(n) fc 



(1)1 




A , 
L -^n,m,l "> j 



(A.15) 



(m) p 



a 



(n) fc 



r a( fc («).?H,9(i)) 



c 



A 



(A.16) 



(m) p 



(0 9 



b 



B Mass Eigenmode Wavefunctions 

Our aim here is to determine the analytic forms of the orthonormal wavefunctions f n (y) and 
g n (y), which are used in (|2.19|) to express the 5D fields A^(x, y) and A^(x, y) in terms of the 
KK mass eigenmodes Af^ (x) and A? n \ 5 (x). In addition, we will derive the transcendental 
equation (|2.24|) that determines the masses for the KK gauge fields. 

We start our discussion with the observation that the gauge-invariant part of the 5D 
QCD Lagrangian ()2.16|) is proportional to [1 + r c 8(y)]. Then, our approach to quantization 
in the gauge can be consistently formulated if this <5(?/)-dependent factor is promoted 
to an overall integral weight of the full quantized action that includes the gauge-fixing and 
Faddeev-Popov ghost terms given in (|2.39jl . In such an approach, it is important to require 
that total derivatives of fields or product of fields with respect to y do not contribute to 
the action, but vanish identically. This condition on the surface terms can be achieved by 
requiring that 

pirR+e fTrR+e 



dy [1 + r c 5(y)] d 5 f n (y) = , / dy [l + r c 5(y)] d 5 g n (y) = . (B.l) 



Bear in mind that our compact space y is defined in the interval (—ttR + e, ttR + e], where 
e is an infinitesimal positive constant which is taken to zero at the end of the calculation. 
Using the convolution methods presented in Appendix C, it is not difficult to show that 
(jB.lj) is sufficient to ensure the vanishing of any total derivative when integrated along 
with the <5(y)-dependent weight [1 + r c 8(y)]. As a consequence of this, the undesirable 
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mixing term (d^A^) (85 A 1 *) in the gauge- kinetic part of the 5D Lagrangian cancels against 
a corresponding term (d^A 5 ) (d^A 11 ) that occurs in the gauge- fixing Lagrangian (|2.39|) . 

In addition to (jB.l|) . the wavefunctions f n (y) and g n (y) have to satisfy the ortho nor- 
mality conditions: 

-irR+e 

dy [1 + r c 5(y)] f k (y)fi(y) = h,i , 

(B.2) 

nR+e 



-ttRA 



dy [l + r c 6(y)]g k (y)gi(y) = 5 kjl 



The analytic form of f n (y) and g n (y) can be fully specified by requiring that the kinetic 
part of the effective 4D Lagrangian for the KK modes A?, and A? n ^ 5 takes on the expected 
form in the conventional gauge: 



+ \ (<V*(n) 5 ) (^) B ) ~\< Wn)5) 2 ■ (B.3) 

The above form of the Lagrangian ()B.3|) is obtained if f n (y) and g n (y) satisfy the simple 
wave equations 

d!fn{y) + mlUv) = 0, (B.4) 
d 2 5 g n (y) + m 2 n9n (y) = 0. (B.5) 

To find the solutions to the above wavefunction equations, we should consistently implement 
the constraints (|B.1J) and (jB.2|) . Notice that our approach is different from the one presented 
in the existing literature [20]. 

Let us now describe how to find the analytic form of the wavefunctions f n {y)- Because 
of the presence of the function 5(y) through the constraints, our starting point is the 
piecewise ansatz 

fn(y) = s ( B - 6 ) 

An sin m n y + Bu cos m n y , for e < y < ttR + e . 



Obviously, our regularized ansatz (|B.6|) is a solution to (jB.4|) . which excludes an infinites- 
imal interval (— e, e) in the neighbourhood of the singular point y = 0. The difference 
between the wavefunctions f n (y) and g n (y) is that the f n (y) is an even function of y, 
whereas g n (y) is an odd one. Moreover, we must assume that both f n (y) and g n (y) are 
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periodic functions of y, i.e. f n (y) = fn(y + 2ttR) and g n (y) = g n (y + 2nR). Within their 
definition interval (— nR+e, irR+e], the periodic condition for the points y = ±ttR amounts 
to the constraints: 

i^-Trtf + e) = f n (irR + e), g n (-nR + e) = g n {itR + e). (B.7) 

Although the above periodicity constraint on f n (y) may be trivial because f n {y) is an even 
function, the one applied to y-odd functions, such as d^f n {y) and g n {y), can provide useful 
information (see discussion below). 

The unknown coefficients A\j\ and Bj ^, as well as the KK masses m n , may be deter- 
mined as follows. From the requirement that f n (y) is an even function, i.e. f n (y) = f n (~y), 
we get the constraint 

A 1 = -A n = A n , B l = B u = B n . (B.8) 
Imposing the first condition of (jB.ljl on the y-odd expression d§f n (y) leads to the constraint: 

d 5 f n (y = 0) = . (B.9) 

To implement (|B.9|) in agreement with the wave equation (jB.4|) for f n (y), we analytically 
define f n (y) in the interval — e < y < e as follows: 

f n (y) = R n cosm n y. (B.10) 

The arbitrary constant R n may be determined by the matching condition 

R n = lim ( ton f n (y)) = lim ( lim Mv)) ■ ( B - n ) 

Notice that the above double limits are well-defined, even though f n (y) has finite jumps at 
the critical points y = ±e. The derivative d 5 f n (y) of the e-regularized function Mv) i s a ^ so 
analytically well-behaved at the singular point y = 0, where S(y) acts, and complies with 
the constraint ()B.9j) . In Fig. El we display the y-profile of the e-regularized function f n =2{y), 
and show its shape in the limit e — >• 0. Here, we should stress again that our regularization 
method consists in taking the limit e — * at the very end of our calculation. 

As we will see below, d$f n {y) is an odd function that shares the same properties 
as g n (y)- We may understand this fact by observing that if the even function f n {y) is a 
solution to the wave equation (|B.4|) . its derivative d$ f n (y), which is an odd function, is 
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-jiR + e -eO e „ 7tR + e 



Figure 3: The y -profile of the e-regularized mass eigenmode wavefunctions f n =2{y) and 
g n= 2{y)- The dotted lines indicate the y-dependence of the wavefunctions in the limit e — > 0, 
in the neighbourhood of the singular point y = [cf. \2.20\) and \2. 21]) J. 

a solution to ()B.4|) as well. Thus, assuming a posteriori that, up to a constant, the odd 
function d§f n (y) equals g n (y), we may impose the second condition of (jB.l|) in the form 

dy[l + r c S(y)} d 5 (d 5 f n (y)) = 0. (B.12) 

-nR+e 

Using the wave equation ()B.4j) for the function f n , we get the new constraint: 

d 5 f n (irR + e) - d 5 f n (-7rR + e) - d 5 f n (e) + d 5 f n (-e) = m 2 n r c f n (0) . (B.13) 
In deriving (|B.13p . we have neglected terms 0(e) that occur in the region —e<y<e. 
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Our ultimate constraint arises from the fact that in addition to f n (y) which is periodic 
by construction, d 5 f n (y) has also to be periodic on the entire y interval, i.e. d 5 f n (y) = 
dsfn(y + 2rrR). Applying the periodicity restriction for the points y = ^ttR + e that lie 
within the definition interval, we get the non-trivial relation 

d 5 f n (-7rR + e) = d 5 f n (irR + e). (B.14) 

The constraints (|B.13j) and (|B. 14j) lead to a determination of the coefficients A n and B n 
and the KK masses m n : 

A n = — - — B n , sinm n 7rit + — - — cosm n nR = . (B.15) 

The second equation in (|B.15|) is equivalent to the transcendental equation (|2.24jl . The 
remaining freedom is an overall normalization constant N n given in (|2.20j) and can be 
determined from the orthonormality condition (jB.2|) . 



To determine the analytic form of g n {y), we proceed very analogously. We start with 
a similar ansatz 

{Cismm n y + D\ cos m n y , for — ttR + e < y < — e 

(B.16) 
C\i sin m n y + _D n cos m n y , for e < y < ttR + e . 

From the fact that g n (y) is an odd function, we get 

Cj = C u = C n , Di = -D u = D n . (B.17) 

In view of the constraints ()B.17|) . g n (y) exhibits a discontinuity in the limit y = ±e — > 0. 
Exactly as we did above for f n (y), we analytically regularize g n (y) in the interval (— e, e) 

as 

9n(y) = Q n sinm n y. (B.18) 

Observe that (|B.18j) is also a solution to the wave equation ()B.5jl for the function g n (y). 
Since dc,g n (y) determined from (jB.16|) has well-defined limits when y = ±e — > 0, we apply 
the matching condition 

m n Q n = lim ( lim d 5 g n (y) ) = lim ( lim d 5 g n (y) J (B.19) 



to determine the constant Q n . Given ()B.18jl . the second equation of the condition (jB.ljl 
implies that 

g n (nR + e) - g n (-nR + e) - g n (e) + g n (—e) = - r c d 5 g n (y = 0) . (B.20) 
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In writing down (|B.20|) . we also made use of the fact that terms 0(e) arising from the 
integration in the interval (— e, e) have been neglected. For illustration, we show in Fig. |3] 
the ^/-dependence of the e-regularized wavefunction g n =2{y), along with its y-profile in the 
limit e —>■ 0. 

Finally, taking into account the periodic condition (jB.7|) for the functions g n {y), we 
obtain the relation 

D n = - C n , (B.21) 

and the second equation in (jB.15|) . As before, the overall normalization constant can be 
determined from the ort honor mality condition (|B.2J) . Finally, implementing the matching 
conditions (jB.ll|) and (jB.19|) . we find that R n = Q n . 

After having derived the analytic forms of f n (y) and g n {y) given in (|2.2U|) and (J2.21j) . 
respectively, it is not difficult to see that f n (y) and g n (y) satisfy the relations in ()2.23j) . 
for the entire piecewise defined region (—ttR + e, ttR + e] that also includes the singular 
point y = 0. Finally, we should remark that for a 5D orbifold theory without BKTs, the 
functions f n (y) and g n (y) go to the standard Fourier expansion modes in limit r c — > 0: 

lim f n (y) = = cos^ , lim . g n {y) = — L= sin ^ . (B.22) 

For completeness, we present in the remainder of this appendix analytic results con- 
cerning the case of a 5D orbifold theory with two BKTs at the orbifold fixed points y = 
and y = ttR. We follow a very analogous approach to the one given above for the one BKT 
case. We start again with the same form of regularized ansaetze f)B.6|) and (|B.16|) for f n {y) 
and g n (y), which are piecewise defined in the interval (—ttR + e, — e] U [e, ttR — e]. Notice 
that the ^/-intervals, (— e, e) and (irR — e, irR+e], are excluded, because of the singular over- 
all integral measure: [ 1 + r c 6(y) + r c 5(y — nR) ] . However, most of the discussion given 
above goes through, with some obvious modifications related to the presence of the second 
5-function, 5(y — ttR). In fact, the only crucial difference is that the constraints (|B.13|) 
and ()B.20|) now become 

d 5 f n (nR-e) - d 5 f n (-7iR + e) - d 5 f n (e) + d 5 f n (-e) = m 2 n r c [f n (0) + f n (nR)} , 

(B.23) 

g n (7rR-e) - g n (-irR + e) - g n (e) + g n (-e) = -r c [d 5 g n (0) + d 5 g n (irR)] . 

(B.24) 

Taking into consideration the above non-trivial modifications, we find that the mass eigen- 
mode wavefunctions f n (y) and g n {y) retain the functional form given in ()2.2()j) and (|2.21j) . 
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in the limit e — > 0, for the one BKT case. However, the normalization factor N n modifies 
in the presence of the second BKT at y = ttR. The analytic form of N n may be determined 
by 



N~ 



1 + 4f c + 



l + 2f c 



+ 



l + 2f e 



1 — r c irRm n ) 2 (1 + r c nRm n ) 2 
1 + 3f c 1 + 3f c 



(B.25) 



1 — r c TiRm n 1 + r c TrRm r 



Moreover, the transcendental equation determining the mass spectrum of the effective 4D 
theory may be cast into the useful factorizable form: 



m. 



tan 



m n TiR 



cot 



m„7rR 



m„r r 







(B.26) 



which is valid for r c ^ and, after some algebra, agrees with [6]. In Fig. we present 
typical solutions to (|2.24|) and (jB.26|) for a 5D orbifold theory with one and two BKTs, 
respectively. 









\ 1-BKT\spectrum • \ 
\ 2-BKT Bpectrum o I 
I i I i 1 



-1/R 



1/R 



2/R 



Figure 4: Mass spectrum m n of the effective J^D theory with one/two BKT in the 5D 
Lagrangian, \2.2Jj\) and \B.26\) . respectively. 
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C Convolution Integrals for S 1 /^ Wavefunctions 

In the process of the KK reduction, expressions that involve products of 5D fields com- 
pactified on an S ,1 /Z 2 orbifold, such as A a Jx,y) and Af(x,y), need to be integrated over 
the compact dimension y. These 5D fields have an even or odd parity under Z 2 transfor- 
mations and can be expressed in terms of the KK mass eigenmodes f n (y) and g n (y), given 
in (|2.2U|) and (|2.21|) . respectively. Using convolution integral techniques, we may express 
the product of two Z 2 -even functions, e.g. F(y) and G(y), in a series of the orthonormal 
wavefunctions f n (y) as follows: 

00 

F(y)G(y) = £ [F * G] {n) f n (y) , (C.l) 

ra=0 

where the coefficients [F * G] r n ) are given by the convolution integrals 

[f * G] (n) = F w G w / d v[ 1 + r Av)] f k (y)fi(y)fn(y) ■ (C.2) 

k,i=o J - jtR 



Likewise, the product of a Z 2 -even function F(y) with a Z 2 -odd one H(y) can be expanded 
in a series of the orthonormal wavefunctions g n {y)'- 

00 

F(y)H(y) = £ [F * H] {n) g n {y) , (C.3) 

ra=0 

where the coefficients [F * H] ( n ) are given by the convolution integrals 

00 pttR 

[F*G} (n) = ^2 F (k)H(i) dy[l+r c 5(y)]f k (y) gi (y)g n (y). (C.4) 
k,i=o J - 7rR 



The integrals ()C.2|) and (jC.4|) appear in the calculation of the effective Feynman rules, 
and of the effective gauge and BRS transformations, when the 5D fields are expressed in 
terms of their KK modes. Specifically, defining the overlap integrals as 



T " A fc . 



dy [l + r c %)]/ fc (y)/,(2/)/n(j/) 



l.n 



-nR 2\ / 2 S ^+ S ^ +S ^ 7r R 
ttR ^ V ' / 

dy [l + r c 5(y)] f k (y) g t (y) f n (y) = , 

we may analytically calculate the coefficients A k t i <n and A^ jn as follows: 

A Khn = A(m k ,m h m n ) + A(-m k ,m h m n ) + A(m k , —m h m n ) + A(m k , m h -m n ) , 

A kin j = - A(m k ,mi,m n ) - A(rn k ,m h -m n ) + A(-m k ,m h m n ) + A{m k , -m h m n ) , 

(C.6) 
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where 



A(m k ,m h m n ) 



N k NiN n 



Jp R dy cos(m fc + + m n )(y - ttR) 
ttR cosm^R cosmi7rR cosm n irR 



+ r c 



(C.7) 



The normalization constants N n are given in (|2.22j) and f c = r c / (2ttR) . According to the 
above definition, the coefficient A k} i >n is symmetric in all of its indices, while A k} i >n is only 
symmetric under the interchange of the first and the last index. Note that the argument 
of the cosine in (|C.7J) may vanish, when some of the arguments in A(m k ,mi,m n ) become 
negative. Taking this possibility into account, we find 



A(m k ,mi,m n ) 



for m k + mi + m n ^ 



N k N t N n { [(1 + 7r 2 R 2 f 2 m 2 )(l + 7r 2 R 2 f 2 m 2 )(l + K 2 R 2 r 2 c m 2 n )\ + f c |, 

for mfc + + m n = . 

(C.8) 

For non- vanishing f c , the lower case is only important for the calculation of A(m k , mi, m n ), 
if at least one of the indices k, I, n is zero and the other two are equal. When this condition 
is satisfied, we find 

4 



Ao,o,o 

Afi,n,0 



y/T+r~ e 



*n,0,n 



2N 2 

n 

y/l + f 







1 + r c + Tr 2 R 2 r 2 ml 



(C.9) 



In any other case (up to symmetries), the coefficients A k ^ n and A k j >n have been calculated 
to have the more compact analytic forms: 



^k,l,n 



^k,n,l 



N k NiN n 7r 2 R 2 fl 



N k NiN n 7r 2 R 2 f 3 



m\ + mf + — 2{m 2 k m 2 + nr[m 2 n + m^ml) ' 
Am k mim 2 (m 2 k + mf — m 2 ) 



(CIO) 



m\ + mf + m\ — 2(mfmf + mfm 2 + m 2 m%.) 
If two of the indices are equal, the last two formulae simplify further to 



A 



n,n,j 



N 2 NjH 2 R 2 rl 



imz 



m A 



Am 2 



A 



Anj^ri 



N 2 N tt 2 R 2 f 3 Am i m n 
n 3 c m 2 -4m 2 n ' 

„ 9 o o 4m5? (m 2 — 2m 2 ) 
— N 2 NjTrRfl 



m 2 — 4m 2 n 



(C.11) 
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Correspondingly, we may now evaluate overlap integrals that involve the product of 
4 orthonormal wavefunctions, 



txR 



dy [l + r c 5(y)] fk(v)fi(v)fMfm(y) 

dy [1 + r c 8(y)] f k {y)fi{y)g n {y)g m {y) 



A 



k,Ln,m 



A k i n m 



(C.12) 



The coefficients &k,i,n,m and A kt i^ m are given by 

^k,i,n,m = A(m k ,m h m n ,m m ) + A(-m k ,m h m n ,m m ) + A(m k , -m h m n , m m ) 

+ A{m k ,m h -m n ,m n ) + A(m k , m h m n , -m m ) + A{-m k , -m h m n , m m ) 
+ A(-m k ,m h -m n ,m m ) + A(-m k , mi, m n , -m m ) , (C.13) 

A fc ,i,n,m = - A(m k ,mi,m n ,m m ) - A(-m k ,mi,m n ,m m ) - A(m k , -mi, m n , m m ) 

- A(-m k , -mi,m n ,m m ) + A(m k ,m h -m n ,m m ) + A(m k ,m h m n ,-m m ) 
+ A{-m k ,m h -m n ,m m ) + A(-m k , m h m n , -m m ) , (C.14) 

where A(m k ,mi,m n ,m m ) is an obvious generalization of (|C.7[) . i.e. 



A(m k ,mi,m n ,m r , 



N k N t N n N m 



2d2~3 mim n m m + m k m n m m + m k niim m + m k mim ri 



7T R T 



m k + mi + m n + m Tt 



x < 



for m k + mi + m n + m m ^ 



+ r ( 



[1 + n 2 R 2 r 2 c m 2 k )(l + n 2 R 2 r 2 c m 2 )(l + 7i 2 R 2 f 2 c m 2 n )(l + n 2 R 2 f 2 c m 2 m 

for m k + mi + m n + m m = . 

(C.15) 



It is easy to see that A ki i^ m is symmetric in all of its indices, while A kt i )n ,m is only symmetric 
in the first pair as well as the last two of its indices. If one of the 4 indices is zero in Ak,i jn ,m 
and A k} i^ m , then (|C.13J) and (|C.14J) reduce to 



A k i n Q 



Ao,k,i 



2A 



l,k,n 



(C.16) 



It is also useful to list analytic results for special combinations of indices, when the lower 
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case in flC.15|) becomes relevant: 



A 



0,0,0,0 



A n n n n 



1 +r c 



2N 2 N 2 

n m 



1 + f c + {l-? c yyR 2 {m 2 n + m 2 m ) + f c Vi? 4 m>™ 



(C.17) 



3A 



n,n,n,n 



f c (l - T^tfrlmi) + (1 + Tr^f 2 ™ 2 ) 2 



Finally, let us comment on the limit of vanishing r c for our analytic results. Taking 
this limit for the expression A(mk,mi,m n ), for example, is a non-trivial task. One should 
take care of the fact that although one may have m k + mi + m n ^ for r c ^ 0, it could be 
that rrik + mi + m n — > for r c — > 0. In such a case, the denominator that appears in the 
first line of (|C.8|) approaches zero, and one has to carefully expand the BKT corrections to 
the KK masses in powers of r c to obtain a sensible result. In this way, we find that 



lim A(m k ,mi,m n ) 

r c — >0 



Analogous considerations for A(mk,mi,m n ,m m ) lead to 



lim A(m k ,m h m n ,m Ti 



k+l+n+m,0 



(C.18) 



(C.19) 



Given that Sk+i+ n ,o and 6k+i+ n +m,o are the usual Kronecker symbols, we have checked that 
our analytic results perfectly agree with those presented in [7], within the context of 5D 
Yang-Mills orbifold theories without BKTs. 



D High Energy Unitarity Sum Rules 

In proving the ET in Section 2, we encounter infinite sums that involve products of the 
gauge-coupling coefficients A fc) ^ n and A k) i ;Tl defined in Appendix C [cf. (jC.6|) ]. For instance, 
we have to carry out the infinite sum, 

oo 
3=0 

that occurs in the calculation of (jH.fijl . Considering (jC.ll|) . the problem reduces to finding 
an expression for the sum, 

^ N 2 

(m 2 — 4m 2 i ) 2 
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The calculation of such infinite sums can be carried out analytically, by extending the 
complex integration analysis techniques developed in [27]. 

Let us briefly describe our complex integration analysis method by considering the 
infinite sum 



oo 
3=1 



where the normalization factors Nj are given in (|2.22|) . We start with the complex analytic 
expression 



z + tamrRz/(r c 7rR) \ -r c cos 2 irRz J z — rrij z — rrij 

where we have expanded it about the regions \z — m,j\ < e, close to its poles ±rrij. Here and 
in the following, we use the convention that m_j = —rrij and rrij > 0, for j > 0. In the last 
equality of (jD.lj) . we have defined the function N(z), in a way such that it is N(m3) = N 2 
at the residuum of (ID. ID. 



Our next step is to use Cauchy's theorem [31] to integrate the LHS of (jD.l|) over the 
poles rrij, i.e. 

lim / dz\z-\ — tan7ri?2;] = lim / dzf c . (D.2) 



The contours C n are the circles z n = (n + 1/2) /R e ie defined in the complex plane, whose 
radii are taken to infinity in a discrete manner, i.e. n — > oo. This ensures that neither 
rrij, nor possible poles of the complex function N(z) do lie on the contour C n . The LHS 
of (jD.2|) then becomes a simple integration over the angle 9 that can be performed easily 
to give 2iri. The RHS of (|D.2|) can be calculated using Cauchy's theorem. Thus, we arrive 
at the equality 

oo 

2m = 2ni f c N{m]) , (D.3) 

j=-oo 

which is equivalent to the desired result 

YN* = ^ . (D.4) 

3 2f c (l + f c ) 1 ; 

The above method can be extended to more complicated infinite sums that include 
factors, such as , : n . For example, a typical infinite sum that arises in such a calculation 
is 

00 N 2 m 4 
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In this case, we may exploit the analytic structure of the summing terms with respect 
to rrij, and consider the analytic continuation rrij — > z in the complex plane. Then, we 
multiply both the LHS and the RHS of (|D.1J) with the analytic expression 

-v4 



(z 2 — 4m 2 ) 2 

This extra factor leaves the integral on the LHS of (jD.2|) intact, since the new factor goes 
to 1 in the limit \z n \ — > oo, for n — > oo. To calculate the RHS of (jD.2|) . we have to properly 
include the contributions from the two additional double poles at z = ±2m n . Thus, we 
finally obtain 

V- N f m j 1 ( 3 \ 

(mj - Ami ? = n**Ii*m£r«N* \ Amn + \ Xn ) ' (D ' 5) 

where A n>nin>n and X n are given in ()C.17|) and ()3.16jh respectively. 

With the help of the complex integration analysis method outlined above, the follow- 
ing list of high energy unitarity sum rules can be derived: 

oo oo 

3=0 3=1 

00 00 

E9A0A A - A \ ^ 9-^,0 A 2 —A -u Y (Ti 7\ 

L - i n,n,j'— l n,3,n L - i n,n,n,n j / ^ n,j,n L -^n,n,n,n 1 ^ v n j v 1 ^' 1 / 

3=0 i=0 



OO 



3=0 3=0 



00 



Eo-5j,o A 2 a , y fP 

n,j,m L - i n,n,m,m ~ 1 n,m j V 1 ^"-'/ 



3=0 

00 

^ ^ 2 3 '° A^jy A n j i?n ^k,l,n,m ^k,l,n,m 1 (P. 10) 

3=0 

where 

K, ro = 4iV n 2 iV^7r 2 J R 2 f c 3 (m2+m 2 n ), (P.ll) 
Z h ,i, n , m = 2N k N l N n N m 7i 2 R 2 r 3 c (m 2 k + m 2 + m 2 n + m 2 J . (P.12) 

The sum rules (IP.6j) and (|P.7jl are relevant to restore unitarity in high energy 2^2 
processes where all KK gauge modes in the initial and final state have the same mass 
m n . The sum rules f)P.8|) and (|P.9|) ensure high energy unitarity for 2 — > 2 scatterings, 
where the KK gauge modes in the initial and final states have the same mass m n and 
m m , respectively, but m n ^ m m . Finally, (jP.10)) is relevant to restore unitarity in 2 — ► 2 
scatterings, when all 4 asymptotic states happen to have different masses. 
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